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Abstract. We consider a natural variation of the concept of stabbing a 
segment by a simple polygon: a segment is stabbed by a simple polygon V 
if at least one of its two endpoints is contained in V ■ A segment set 5 
is stabbed by V if every segment of S is stabbed by V ■ We show that 
if S is a set of pairwise disjoint segments, the problem of computing 
the minimum perimeter polygon stabbing S can be solved in polynomial 
time. We also prove that for general segments the problem is NP-hard. 
Further, an adaptation of our polynomial-time algorithm solves an open 
problem posed by LofHer and van Kreveld [Algorithmica 56(2), 236-269 
(2010)] about finding a maximum perimeter convex hull for a set of 
imprecise points modeled as line segments. 



1 Introduction 

Let <S* be a set of n straight line segments (segments for short) in the plane. 
The problem of stabbing S with different types of stabbers (in the computer 
science literature) or transversals (in the mathematics literature) has been widely 
studied during the last two decades. 

Rappaport |16j considered the case in which the stabber is a simple polygon. 
Specifically, he studied the following problem: a simple polygon V is a polygon 
transversal of S, if we have V l~l s ^ for all s 6 S; that is, every segment 
in S has at least one point in V. A simple polygon V is a minimum polygon 
transversal of S if V is a polygon transversal of S and all other transversal 
polygons have equal or larger perimeter. Rappaport observed that such a polygon 
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always exists, it is convex, and may not be unique. He gave an 0(3 m n + nlogn)- 
timc algorithm for computing one, where m is the number of different segment 
directions. Several approximation algorithms are known |5I8| . but determining if 
the general problem can be solved in polynomial time is still an intriguing open 
problem. 

Arkin et al. [2] considered a similar problem: S is stabbable if there exists 
a convex polygon whose boundary C intersects every segment in S; the closed 
convex chain C is then called a (convex) transversal or stabber of S. Note that in 
this variation there is not always a solution. Arkin et al. [2] proved that deciding 
whether S is stabbable is NP-hard. 

In this paper we also consider the problem of stabbing the set S by a simple 
polygon, but with a different criterion that is between the two criteria above. 
More concretely, we use the following definition: 

Definition 1. A segment s G S is stabbed by a simple polygon V if at least one 
of the two endpoints of s is contained in V. The set S is stabbed by V if every 
segment of S is stabbed by V . 

With this definition we study the Minimum Perimeter Stabbing Poly- 
gon (MPSP), defined as finding a simple polygon V of minimum perimeter that 
stabs a given set S of segments. The MPSP problem is radically different from 
the two problems above, those studied by Rappaport [TB] and Arkin et al. [2], 
because for the MPSP only the endpoints of the segments play a role in the 
solution. Indeed, an alternative way to describe the input to the MPSP problem 
is by saying that the input are pairs of points instead of segments. However, as 
we will show in this paper, the segments play an important role in establishing 
the difficulty of the problem, hence we stick to the original definition. 

Moreover, the difference with the problem of Rappaport [TB] is that in his 
definition V can have both endpoints of a segment of s € S not in V (provided 
that the interior of s is stabbed by V), whereas we force one of the endpoints 
to be in P. One of the common properties of both problems is that the optimal 
solution is a convex polygon and that it always exists (the convex hull of S is 
always a stabbing polygon). 

On the other hand, a difference with the definition used by Arkin et al. is 
that in the MPSP problem a segment of S can be fully contained in V, with 
both endpoints in the interior of V, while this is not allowed in the problem 
studied by Arkin et al. Therefore, we can say that our problem is between the 
two mentioned ones. 

Related work. Prior to the paper by Rappaport [HI], Meijer and Rappaport [T5] 
solved the same problem for a set of n parallel segments in optimal <9(rtlogn) 
time. Mukhopadhyay et al. |15j considered a similar problem in which the seg- 
ments are all vertical, and proposed an 0(n log n) time algorithm to find a 
minimum-area convex polygon transversal of S. For parallel segments, Goodrich 
and Snoeyink [7| gave an 0(n log n)-time algorithm that decides whether a con- 
vex transversal exists. 

Several similar problems have been considered in the context of data impreci- 
sion by Loffler and van Kreveld |12|13j . Their input is a set of imprecise points, 
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where each point is specified by a region in which the point may lie. The output 
is the smallest and the largest possible convex hulls, measured by perimeter and 
by area. Among the results obtained in |12j , we cite those where regions are seg- 
ments. For maximum-area convex hulls, the problem can be solved in 0(n 3 ) time 
if the segments are parallel, or when they are pairwise disjoint with endpoints 
in convex position. The problem is NP-hard for general segments. 

The minimum-perimeter and minimum-area convex hulls problems for par- 
allel segments coincide with the problems studied by Meijer and Rappaport [2] 
and Mukhopadhyay et al. [15] . respectively. Notice also that the setting we con- 
sider is in fact a constrained version of the problems studied by Loffier and van 
Kreveld |12j . in which each imprecise point is specified by a pair of points. 

Pairs of points are also the input to the problems studied by Arkin et al. pQ , 
who studied the 1-center and 2-center problems for pairs of points. In the former 
problem, the goal is to find a disk of smallest radius containing at least one point 
from each pair. The latter one aims at finding two disks of smallest size such that 
each pair has one point in each disk. Arkin et al. [I] presented algorithms for 
these problems that run in O(n 2 polylog n) and 0(n 3 log 2 n) time, respectively. 

In a more general setting, Daescu et al. [I] studied the complexity of the 
problem of given a fc-colored point set, finding a convex polygon of minimum 
perimeter containing at least one point from each color. Note that the MPSP 
problem is the special case in which 2n points are colored with n colors and each 
color is used twice. They proved that their problem is NP-hard if k is part of 
the input of the problem. 

Our results. We show in Section [2] that if S is a set of pairwise disjoint seg- 
ments, the MPSP problem for S can be solved in polynomial time. We then 
show how the algorithm can be adapted to solve the following maximization 
problem: Select exactly one point on each segment in S such that the perimeter 
(or area) of the convex hull of the selected points is maximized. This problem 
was stated as open |12) . and is also the solution to the maximization variant of 
the transversal problem [12]. In Section [3] we show that for general segments the 
MPSP problem is NP-hard. We complement the NP-hardness by showing that 
the MPSP problem is Fixed Parameter Tractable (FPT). 

Note throughout the paper that optimization on the perimeter requires com- 
paring sums of radicals (specifically, the sum of Euclidean distances). It is not 
known whether this problem is in NP [3] , and therefore the NP-hardness result 
does not imply NP-completeness for the decision version of the problem. For the 
same reason, we assume the real RAM as the underlying computational model in 
our algorithms. Since our algorithms are combinatorial and only the cost func- 
tion depends on the geometry of the problem instance, the methods in Section [2] 
are also applicable for optimizing the area (which is in NP). 

2 Solving the problem for pairwise disjoint segments 

In this section we show that if the segments in S are pairwise disjoint, then the 
MPSP problem can be solved in polynomial time. Given any two points p and q 



3 



in the plane, let pq denote the segment joining p and q. For any simple polygon V 
let dV denote the boundary of V . 

Consider all possible bitangents of S, i.e., let B be the set of all segments 
not contained in S spanned by two endpoints of segments in S. Note that the 
elements of B might cross each other and might also cross the segments in S. 
A polygon C* with minimum perimeter that contains at least one endpoint of 
every segment of S is spanned by endpoints of segments in S, and its edges are 
elements of B. 

Arkin et al. [2] describe a dynamic programming approach to decide whether 
a set of pairwise disjoint segments admits a convex transversal (the vertices of the 
transversing polygon are restricted to a given set of candidate points) . They use 
constant-size polygonal chains that separate subproblems and are not crossed by 
segments; therefore the subproblems are independent. We adapt their approach 
to produce an algorithm for the MPSP problem. The main difference (apart 
from the fact that no candidate points are needed) is that segments actually can 
cross the separating chains. However, we show below that they can be handled 
in a way that leads to polynomial running time. Afterwards, we discuss how to 
adapt this approach for the maximization variation. 

Triangulating a combination of segments and a polygon. The following 
way of triangulating a combination of segments and a polygon is crucial for the 
algorithm, and motivates the structure of the subproblems used in the dynamic 
programming algorithm. 

Let Q be a simple polygon and let S c be a set of pairwise disjoint segments 
of which each crosses dQ exactly once. Note that throughout this section we 
distinguish between a segment intersecting (having a point in common) and 
crossing (having an interior point in common with) another segment or set. 
Let X be the interior of Q and let X' denote the set we get after removing the 
1-dimensional domains of S c from X, i.e., X' = X\\j seS s. Then X' is an open 
region whose closure is Q. Note that the vertices of X' are the union of: (i) the 
vertices of Q, (ii) the endpoints of edges in S c that are in the interior of Q, and 
(iii) the points where elements of S c cross dQ. Further, note that X' might not 
be connected if there is a segment of S c that has one endpoint on dQ and the 
other one outside Q (e.g., the longest segment in Fig. [I] left). 

We now triangulate X' (i.e., partition it into triangles that are spanned 
only by vertices of X' , see Fig. [I]). The triangulation T of X' behaves like the 
triangulation of a collection of simple polygons (imagine the 1-dimcnsional parts 
not in X' where the segments of S c enter Q, i.e., X\X' , to be slightly "split", as 
in Fig.JT] center). Note that the vertices of T are exactly the vertices of X' . Each 
edge in T that is not part of dQ or part of a segment in S c partitions X' into 
two sets (note that each set need not be connected). We call such edges chords 
(gray edges in Fig. [TJ right). Chords are the equivalent of diagonals of simple 
polygons (interior edges that subdivide the polygon into two smaller polygons). 
Further, X' might also be separated by an edge that is part of a segment in S c 
(like the longest edge in Fig. [T]) . We call such a segment a separating segment. 
Keep in mind that there are chords that have one or both of their endpoints 
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Fig. 1. Left: an optimal polygon Q, only the solid edges are in S c - Center: schematic 
view of X' as a collection of simple polygons. Right: a triangulation of X' , gray edges 
are chords. The segments fully contained in the polygon (shown dashed) are ignored 
by the triangulation. 




Fig. 2. Examples of bridges. The two bitangents defining the subproblem are shown 
dashed, chords are dash-dotted, and segments from S c are shown solid. 



not on the endpoint of a segment or at a vertex of Q, but at the crossing of 
a segment with dQ. In any case, a chord or a separating segment defines a 
polygonal path from one point on an edge of Q to another point on an edge 
of Q. Following [2J, we will use these polygonal paths of at most three edges, 
called bridges, to define our subproblems to obtain a solution when taking the 
MPSP C* as Q. One may think of the approach being similar to the classic 
dynamic programming algorithm for minimum weight triangulations of simple 
polygons |6lllj . but with a major difference: we do not know the boundary of 
the triangulated region beforehand. 

Subproblems. Every subproblem is defined by an ordered pair (a, b) of directed 
bitangents of B and a polygonal chain (3 of at most three edges, the bridge, which 
connects a and b. When evaluating a subproblem (a, b, j3), we assume that a and b 
are edges of C* and that C* equals Q in the discussion above (for some choice 
of S c to be defined later). Therefore, the bridge j3 is part of a triangulation of X' 
and separates X'; (3 is either a part of a separating segment or consists of a 
chord (called the chord of j3) and at most two parts of segments of S c . See Fig. [2] 
for examples of bridges. Note that a bridge might have a chord that is not a 
bitangcnt of B (like the second from the left in Fig. [2]). Further, note that a 
bridge can only be crossed by a segment through the chord, since the segments 
are pairwise disjoint by definition. 

Let the directed bitangents be a — a\02 and b = 6162- Given a directed 
bitangent a — a\ai we write a for the directed bitangent 020,1. W.l.o.g. let a\ 
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Fig. 3. Examples of subproblems. Rightmost: example for the initial pair. 



and bi be on the x-axis and ai and 62 be above it. Also, let b be to the left of 
the directed line through ai and et2. See Fig. [3] for an illustration. 
Solution of a subproblem. We define the solution of a subproblem as follows. 
Let C* jobea polygon of minimum perimeter that: (i) contains a and b as two 
of its boundary edges, (ii) contains at least one endpoint of each segment in S, 
and (iii) contains both endpoints of every segment of S that properly crosses (3. 
The third condition is particularly important, as will become clear later. 

Let C a ,b,i3 be the polygonal chain on dC* b ^ starting at ax, counterclockwise 
traversing dC* b p and ending at b±. Note that C a ,f>,/3 is an open polygonal chain, 
as opposed to C* b/3 , which is a simple polygon. 

The solution of a subproblem (a, b,f3) is C ay b./3, and its cost is the length of 
that chain. The base case occurs when a 2 ~ 62, and has cost equal to the sum 
of the lengths of a and b. Note throughout the construction that this is the only 
way a and b can intersect. In general, a and b form a quadrilateral a,2a\b\bi. If 
the quadrilateral is not convex, we discard the subproblem (i.e., we assign it a 
cost of 00). The general case where it is convex is discussed next. 
Outline of the algorithm. From now on we assume that a and b define a 
convex quadrilateral. The outline of the algorithm is as follows. We guess a pair 
x,y S B such that yiy\X\Xi are four consecutive vertices of C* . Hence, after 
OdS"! 4 ) guesses we have found x and y such that dC* = C XtVt p U y\X\ with 
A) = X \V\ ■ Suppose we are given the solution Q — C* . Let X' be defined as 
above, and let S c be the set of segments in S that cross C x ^ y ^ (which does not 
include the ones that cross /3q). Let A$ be the triangle of a triangulation T of 
X' that has fio — y\X\ as one side. The subproblem (x,y,0o) will be solved by 
guessing the third endpoint of Aq and the edge c of C x _ v ^ that is incident to 
Aq or that is crossed by a segment whose endpoint is incident to Aq. In the most 
general case, this gives two new subproblems (x, c, j3i) and (c, y, A), where each 
of Pi and /?2 contains one side of Aq that is not part of /3q (we will consider the 
other cases in detail below). See Fig. |3j right. 

Let a be the ray through 0,2 starting at a\. Let b be defined analogously. For 
every subproblem (a, 6, /3), only a part of the elements of S is relevant. Consider 
the (possibly unbounded) maximal region to the left of a and to the right of b 
(recall that a and b are directed). The bridge (3 disconnects that region into two 
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parts. The subproblem region R a ,b,f3 is the part "above" /3 (i.e., the part adjacent 
to a \ a and b \ 6; the bridge might not be x-monotone) . 

The subproblem region is marked gray in the examples of Fig. [3] Naturally, 
only the segments that have at least one endpoint in R a ,b,/3 are relevant for 
finding C a ,b,/3- We distinguish between three different types of such segments: (1) 
Segments that are entirely inside R a ,b,p are complete. (2) Segments that share 
more than one point with R a ,b,p but are not complete are cut. (3) A segment 
with infinitely many points on the bridge is neither cut nor complete. We say 
that a point is inside C a ,b,i3 when it is contained in the closure of the region 
bounded by C a ,b,/3 and (3. 

If there is a segment that is entirely to the right of a or to the left of b, 
then the choice of a and b cannot give a solution and such a subproblem is 
assigned oo as cost. We also do this if a segment intersected by a or b does not 
have an endpoint inside the subproblem region. 

Observe that if a segment in a valid subproblem intersects a or b, then we 
know which of its endpoints must be inside C a> b,/3, while we do not know that 
for the cut segments that intersect the chord of the bridge. However, we will 
choose our subproblems in a way such that all endpoints of cut segments in the 
subproblem region will be inside C 0j b iy g; the reason for that will become clear in 
the proof of Lemma [3j but the reader should keep this in mind as an essential 
part of the method. For the complete segments, we need to decide which endpoint 
to select. 

Lemma 1. Given a subproblem instance (o, 6, /3), let t be the chord of P, or its 
only edge if ft is a single edge ( which may be a chord itself, or part of a separating 
segment). Let X be the region bounded by C a ,b,p U (3, and let X' = X \ [J s€ g s, 
for S c the set of segments of S that are crossed by chain C a ,b,p- Then either t is 
an edge of C a ^^, or there exists a triangle A such that: 

1. The interior of A is completely contained in X' . 

2. The edge t is an edge of A. 

3. The apex of A (i.e., the vertex not on t) is either (i) an endpoint of a segment 
in S c inside X , (ii) an endpoint of a segment in S that is a vertex of C a ,b,i3, 
or (Hi) an intersection point between a segment in S c and C a .b,p- 

Proof. Arbitrarily triangulate X' .1ft is not on the boundary, then the triangle A 
incident to t inside the subproblem region fulfills the properties. See Fig. [4] □ 

Lemma 2. Let A be the triangle of Lemma [7} Any segment of S that has a 
non-empty intersection with the interior of A either has both its endpoints inside 
C a ,b,p or crosses t; in the latter case the endpoint that is inside R a ,b,p is also 
inside C a ,b,p- 

Proof. This follows from the properties of A in Lemma[TJ A segment intersecting 
the interior of A is not part of S c but has a non-empty intersection with X. 
Therefore, either both of its endpoints are inside C a ,b,fi, or it enters X via t and 
therefore has its relevant endpoint inside C a ,b.fi by definition. See Fig. [4] □ 
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Fig. 4. Illustration of Lemma [T] Left: four possibilities for A shown in gray. C a ,b,0 is 
dash-dotted, with the defining bitangents dashed. Right: a triangulation of X' . 

Getting smaller subproblems. Let A be the set of points that are either 
endpoints of S or crossing points of a segment and a bitangent (recall that 
no segment of S is an element of B). Hence, A contains all the points that 
are possible apices for a triangle A of Lemma [l] Note that one may construct 
subproblems where every possible apex of A is an endpoint of a segment in S c , 
as well as subproblems where every possible apex is on a point where a segment 
crosses C a>ht p. Further, note that |A| e 0(\S\ 3 ) since \B\ = 4('f l ). 

Consider again the subproblem (a, b, fa}. As in Lemma[TJ let t be the chord of fi 
if a chord exists, or let t otherwise be the only edge of /3. Let ap be the intersection 
point of a with the bridge fa, bp is defined analogously. For each subproblem 
(a, b, fa that is not a base case (i.e., a 2 ^ b 2 ), one of the following cases applies, 
allowing to get one or two smaller subproblems. During the execution of the 
algorithm we will consider both cases. 

Case 1: t is an edge of the solution, i.e., an edge of C a ^ t p. This happens 
when t is a chord that does not intersect the interior of the quadrilateral defined 
by a and b. This case is only valid if no segment crosses t, as we require all the 
endpoints in R a ,b,/3 of segments crossing t to be inside C a ,b,p- In that case we get 
at most two new subproblems (a,t, fa) and (t, b, fa), where fa is the edge apt\ 
and fa is the edge t 2 bp. However, note that one of (a,t) or (t, b) (or both) might 
intersect at a 2 or b 2 , respectively, and therefore form a base case. 

Case 2: t is not an edge of the solution. Then there is a triangle adjacent 
to t as in Lemma [l] We will guess the apex of the triangle. For every point d in 
A n R a ,b,0 consider the triangle Ad that d forms with t. We only consider d if 
Ad is completely inside R a .b,p, and where the interior of Ad does not intersect 
any segment that intersects a or b. It follows from Lemma [T] that one of triangles 
tested leads to a subdivision of the optimal solution. We get the following two 
subcases, see Fig. [5] 

Case 2.1: d is a point where a bitangent and a segment cross. Let c be 
the bitangent that contains d. If c equals a or b, then we get one new subproblem 
(a, b, fa), with fa containing a side of Ad as a chord (Fig. [5]i). Otherwise, we get 
two new subproblems, (a,c,fa) and [c,b,fa), where fa and fa both contain a 
side of A d (Fig.§)). 

Case 2.2: d is an endpoint of a segment. Let s be the segment that has 
d as its endpoint. Choose a point x where s intersects some bitangent c. Then, 
for every possible choice of x (which implies the choices of c), we get two new 
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Fig. 5. Case 2. The new bridges are dotted, (a)-(b) Case 2.1. (c)-(d) Case 2.2. 

subproblems (a,c, /3i) and (c, 6, ,$2), as in the previous case; note that for both 
new bridges, x = d is possible. The degenerate case where c equals a or 6 can be 
handled as in the previous case. See Fig. [5j>d. 

Lemma 3. Given any valid subproblem (a,b,/3), there is a pair of subproblems 
among the ones above such that the union of their solutions is equal to C a ^.j3- 

Proof. Consider the edge t of LemmaJT] If t is a chord and part of C a .b,p, then it 
will be considered in Case 1. Otherwise, consider the triangle A inside C a j,p- All 
segments that are intersected by the interior of A are either completely contained 
in C a ,b,/3 or enter through t (if it is a chord) and therefore have their relevant 
endpoint inside C a ,b.p (cf. Lemma [2]). Hence, when the choice of Ad coincides 
with A, the two subproblems can be combined to C a ,b.p\ the only segments that 
are part of both subproblems intersect the interior of A, and we know that both 
endpoints will have to be inside the chain that results from the combination of 
the solutions of the subproblems. Since all possibilities of Ad are checked, the 
subproblem combination of minimum cost is guaranteed to be C a .b^- ^ 

This last lemma now implies that we actually find the optimal solution. Note 
that it is easy to construct a pair of bitangents and a bridge (a, b, /3) that is 
part of the optimal solution but for which C Ui b t p is not part of C* . However, as 
mentioned in the outline of the algorithm, we choose the initial problem (x, y, /3q) 
in a way that dC* — C x>Vi p U j3 . All segments crossing (3 — X\y\ need to have 
their endpoint above /?o inside the solution, and the algorithm actually produces 
a triangulation of X' when taking C* as Q and S c being the segments that cross 
dC* but do not cross /3 . 

Recall that we initialize the algorithm using a brute-force approach: that 
is, we consider all Od^ 4 ) possible choices for two defining bitangents and a 
bridge a\bi. Every subproblem contains less edges of the complete graph on all 
endpoints of S, and for every subproblem we need polynomial time. The number 
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of subproblems can be bounded by the choices for c and d. Therefore, dynamic 
programming can be applied to obtain a polynomial-time algorithmF] 



Theorem 1. Given a set S of pairwise disjoint segments, a Minimum Perime- 
ter Stabbing Polygon (MPSP) — i.e. a minimum perimeter polygon containing at 
least one endpoint of each segment in S — can be computed in polynomial time. 

Maximization for pairwise disjoint segments. While in the problem setting 
of Arkin et al. [2] the boundary of the polygon has to intersect all segments, 
our problem setting only requires at least one endpoint of each segment to be 
contained in the solution. Our previous algorithm relies on the fact that the 
result has minimum perimeter: this automatically prevents two endpoints of the 
same segment from being vertices of the resulting polygon. However, making 
the algorithm slightly more sophisticated, we can solve in polynomial time a 
maximization version of the problem: select exactly one point on each segment 
in S such that the perimeter (or area) of the convex hull of the selected points 
is maximized; stated open by LofHer and van Kreveld [12] . This result is based 
in the fact that for the maximum area or perimeter transversal, one needs to 
consider only the endpoints of the segments jl2, Lemmata 1 and 8]. 

Theorem 2. There exists a polynomial-time algorithm that selects exactly one 
point on each segment in S such that the perimeter (or area) of the convex hull 
of the selected points is maximized over all possible selections. 

Proof. Due to [HI Lemmata 1 and 8], we know that we only need to consider the 
endpoints of the segments. We modify the algorithm used for the minimization 
version of the problem. Note that the structure of the solution is very similar. 
Again, let C* be the optimal solution. One main difference is that a segment 
that has an endpoint as a vertex of C* might have the other endpoint in the 
interior of C*, i.e., might be completely contained in it. We define subproblems 
and bridges in the same way. The crucial property in the previous algorithm was 
that a segment that enters a subproblem region through the chord of the bridge 
had to have its endpoint that was inside the subproblem region to be inside the 
solution of the subproblem as well, motivated by Lemma [2] and the choice of the 
initial bridge (3q. We did not distinguish between the two types of segments that 
enter the subproblem region through the chord, but now we need to do so. 

There are (i) the ones that cross /3o\ we denote these by Sq. Further, there 
are the ones that (ii) are completely contained in the interior of C*; let the set 
of these segments be called S- m t- Recall the proof of Lemma [I] If we replace S c 
by the set of the segments that intersect C a ,b.j3 but are not in Sq, the analogous 
result follows. 

Following the proof of Lemma [2j we observe that the segments of S^t have 
both endpoints in the interior of the solution, while a segment in Sq might have 

5 A straightforward analysis of the running time of the algorithms gives Od^ 9 ), which 
probably can be improved. In any case, it is worth stressing that our main contribu- 
tion is that the problem can be solved in polynomial time, more than the running 
time itself. 
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the corresponding endpoint on the boundary of the subproblem solution. There- 
fore, the choice of the bitangent c that gives new subproblems for a subproblem 
(a, 6, /3) can be altered in the following way. 

If c shares an endpoint with a segment that has its other endpoint on a or 6, 
then c is not valid. Further, c must not share an endpoint with a segment that 
crosses /3 and is not in Sq (however, the requirement that all endpoints in R a ,b,/3 
of segments that cross (3 have to be inside the subproblem solution persists) . 

The last modification is that in Case 2.2, x can also be the endpoint of 
the segment that is not d (recall that the solution might completely contain a 
segment that contributes a vertex to it). 

With this variation, we never select both endpoints of a segment but still 
find (a triangulation of) the optimal solution. □ 

3 Hardness of the general version 

In this section we prove that the MPSP problem is NP-hard by reducing 3-SAT 
to it. Our reduction is similar to the ones used in |2l4ll2j . 

Theorem 3. The MPSP problem is NP-hard. 



Fig. 6. (a) Overview of the reduction from 3-SAT. Variable gadgets (b) are to the left 
and clause gadgets (c) to the right. 

Proof. Let a 3-SAT instance consist of n variables x\, . . . ,x n and m clauses 
Ci, . . . , C m . We reduce this instance to the following one of the MPSP problem. 
We draw a circle and place variable gadgets in the left semicircle, clause gadgets 
in the right semicircle, and segment connectors joining variable gadgets with 
clause gadgets. See Fig. [6^,. 

Gadgets. For each variable Xi, i £ [l-.n], we put points Ti and Fi on the 
circle and place three segments: segment Tji 7 ^, and two zero-length segments Oj 
and bi, so that TiFi is parallel to the line containing both a, and hi. Refer to 
Fig. [7]d. Furthermore, trapezoids with vertices cii,Ti,Fi,bi, for all i £ [l..n], are 



congruent. Let P v := \a t T t \ + |T<6<| = \cnFi\ + and P' v := \ ai T,\ + |T 4 ^| + 

\Fibi\. 
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Fig. 7. Repetition of Fig. [6] (a) Overview of the reduction from the MAX-E3-SAT 
problem. Variable gadgets (b) are to the left and clause gadgets (c) to the right. 

For each clause Cj, j € [l..m], we first place two zero-length segments Cj and 
dj. We select the three points Pji, Pj,2, and Pj,3, dividing evenly the smallest 
arc of the circle joining Cj and dj into four arcs, and then we place three other 
segments: Pj,iPj,2, Pj,2Pj,3, and Pj,3Pj,i- See Fig. [Tj;. The convex pentagons with 
vertices Cj, dj,Pj,i,Pj t 2,Pj,3, for all j G [l..m], are congruent. Let P c :— \cjPj t i\ + 
\Pj,iPj,2\ + \Pja d j\ = \cjPj,i\ + \Pj,iPj,n\+\Pj,a d j\ = \ c jPjM + \Pj,2PjA+\Pj,3dj\ and 
P' c := \cjPj t i | + |Pj,iPj,2 I + \Pj,2Pj,3 1 + \Pj,3<lj | ■ We further ensure that m(P' c —P c )< 
P' u — P v . This condition will be necessary in the problem reduction. 

For each clause Cj, j € [l..m], we add segments called connectors as follows. 
Let Xi be the variable involved in the first literal of Cj . If appears in positive 
form then we add the segment TiPj t \. Otherwise the segment FiPj^i is added. We 
proceed analogously with the variable in the second literal and point Pj,2, and 
with the variable in the third literal and point Pj,3- 

Problem reduction. Consider the set of segments added at variable gadgets, 
clause gadgets, and connectors as an instance of the MPSP problem. Observe 
that any optimal polygon V op t for this instance satisfies the following conditions: 

(a) V op t contains as vertices points et^ and for all variables Xj, i € [l-.n], and 
points Cj and dj for all clauses Cj, j € [l..m]. 

(b) For each variable Xj, i G [1 •■"■], Vopt contains exactly one of Tj and Fi as 
vertex between aj and 6j. 

(c) In the clause gadget of each clause Cj, j € [l..m], if the selected endpoint of 
at least one connector is not in the gadget as a vertex of Vopt, then exactly 
two points among pj t ±, Pj t 2, and Pj t 3 are vertices of V op t- Otherwise, all three 
are vertices of V opt - 

Condition (a) is trivial, condition (b) is due to m{P' c — P c ) < P' v — P v , and 
condition (c) is a consequence of segments Pj,iPj,2, Pj,2Pj,3i and Pj,3Pj,i- Refer 
to Fig. [8] for condition (c). 

Let V be any feasible polygon satisfying conditions (a)-(c). Polygon V induces 
the following assignment for variables x\,X2, ■ ■ ■ ,x n : we assign x, to true if point 
Ti is a vertex of V, false otherwise. With this assignment we can ensure that every 
clause Cj is satisfied if and only if exactly two points among pj t i, Pj,2-, and pj t 3 
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Fig. 8. In each clause gadget, if at least one connector has an endpoint not in the 
gadget as a vertex of optimal polygon Vopt, then exactly two points among Pj,i, Pj,2, 
and are vertices of V op t (any of the first three figures from left to right). Otherwise, 
all three are vertices of V op t (rightmost figure). 



are vertices of V . Therefore, the formula is satisfiable if and only if the perimeter 
of the MPSP is |&id| + \a n d m \ + \a t b l+1 \ + Y^=i \djC j+1 \ + nP v + mP c ). 

The coordinates specifying positions of the points of the gadgets are shown to 
be rational and polynomial in n and m in Appendix [A] Note that our reduction 
uses segments of zero length. This can be changed by locating an endpoint 
sufficiently far from the circle where the gadgets are placed on. □ 

Observe that the same reduction with minor modifications applies for the 
case of minimizing the area of the output polygon. Moreover, our proof shows 
that the problem remains NP-hard even if the endpoints of all the segments 
are in convex position. On the other hand, the -s/2-approximation algorithm of 
Daescu et al. [3] gives the same approximation ratio for our MPSP problem. 

It is worth mentioning that the MPSP problem is FPT on the number k 
of segments that intersect other segments. Namely, let S' C S be the set of 
segments of S that do not intersect with any other segment of S, and consider 
all the 2 k instances of the MPSP problem such that the input of each of them 
consists of the elements of S' joint with exactly one endpoint (i.e., a segment 
of length zero) of each element of S \ S' . All these instances can be solved in 
0(2 k P(n)) time since all of them consist of pairwise disjoint segments, where 
P(n) is the polynomial time of Theorem [l] The optimal solution of the MPSP 
problem for S is among the 0(2 k ) solutions found for those instances. 
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A Exact point construction for Theorem [3] 

In this section we give the exact construction for the segment set described in the 
proof of Theorem [3] We place the segment endpoints of the gadgets at rational 
coordinates on the unit circle in such a way that the size of both the numerator 
and denominator of each coordinate arc bounded by a polynomial function of 
the size of the initial problem. 

With respect to the segment endpoints, both the variable gadgets and the 
clause gadgets have the same structure. For any point p, let Zp denote the polar 
angle of the point. If we construct a variable gadget such that a v := Za; — ZTj = 
ZTi — ZFi = ZFi — Zbi, the gadget works as described. If we use the same angle 
a v for all variable gadgets, the gadgets are congruent. The analogous holds if we 
choose an angle a c — Zcj — Zpji = . . . = Zpj^ — Zdj between two consecutive 
points when constructing a clause gadget; we obtain a set of congruent gadgets 
that fulfill the properties described previously. In the remainder of this section 
we show how to choose the reference points a* and Cj for each gadget among the 
rational points on the unit circle and the angles a v and a c . 

We use several well-known facts about rational points on the unit circle (see, 

e-g-j US])- For any t € Q, the point p t = (^jt^j, t^pij i s rational and lies on the 
unit circle. Hence, the coordinates of pt describe the cosine and sine, respectively, 
of the polar angle Zp t . Observe that for any t > 1, point pt lies in the first 
quadrant]^] In particular, we have Zp t — > and pt — > (1, 0) when t — > 00. Using 
the trigonometric identity sin(ai — a%) — sin(ai) cos(a2) — cos(ai) sin(o!2) we 
obtain: 

Mz Pt - z Pt+l) = (friFTl) 

_ 2(t 2 + t + l) 
~ (i 2 + l)(t 2 + 2i + 2) ■ 

Observation 1 For t > 1, the angle Zp t — Zp t+ i is monotonically decreasing 
in t. 

The following inequality will allow us to choose both the reference points and 
the small angles for the gadget construction. If 1 < t < 5N for some positive 
integer TV (the factor 5 is chosen with foresight), we have 



t 2 -l \ ( 2(t + l) \ 



2(257V 2 + 57V + 1) 

sm(Zpt - Zpt+i) > 



> 



(257V 2 + l)(25iV 2 + ION + 2) 
5(W 2 5(W 2 



(25iV 2 + l)(25iV 2 + 10A^ + 2) 6257V 4 + 25(W 3 + 75iV 2 + 10A^ + 2 

5(W 2 50 



> 



6257V 4 + 250iV 4 + 75A^ 4 + 10iV 4 + 27V 4 962iV 2 



Note that there are multiple conventions for choosing the sign of the coordinates in 
this parametrization. 
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We can use this bound to define a rational angle Zp s that is smaller than all 
intervals we consider in the construction: 



sin(Zp s ) 



2s 




50 



,s 2 + l 



< 



962iV 2 



which is fulfilled if 



s > 



9627V 2 
25 



Let us place the endpoints for the variable gadgets. We choose a, = P(5i_4) 
and therefore set N = n. Further, we choose s = 100n 2 , which fulfills the above 
inequalities. By the choice of s, we have Za i+1 + 3Zp s < Zai, hence the variable 
gadgets do not interfere with each other. We therefore can place T i} Fi, and 
bi on the arc segment between ctj and a,+i by rotating aj up to three times 
by a v = Zp s . The points can be explicitly computed from dj by using the 
coordinates of p s as elements of the rotation matrix: 



for k € {1,2,3}. Observe that the coordinates of a, and p s are the sines and 
cosines of the corresponding angles and are rational; therefore, the resulting 
points also have rational coordinates, which are bounded by a polynomial func- 
tion in the size of the input, since k is constant. Finally, we change the signs of 
the coordinates, so that the variable gadgets are placed on the third quadrant. 

For the clause gadgets, we can basically proceed in the same manner, choosing 
N = m in the above equation, as well as Cj = P(5j_4). However, there is an 
additional constraint to fulfill: m{P' c — P c ) < P' v — P v . This is again an inequality 
in the same style as for the variables and thus also results in rational coordinates 
polynomial in n and m with a parameter s' . Replacing N by m in the previous 
calculation we get the constraint s' > 96 ^™ for the angle a c = Zp s > . The second 
constraint m{P' c — P c ) < P' u — P v is more sophisticated. On cither side of the 
inequality we subtract two Euclidean distances, in particular twice the distance 
between two neighboring points minus the distance between the two neighbors 
of a point. This is given by 



where t = s in P' v — P v and t = s' in P' c — P c . Hence, we need to assure that 





m 



s' 2 + 1 



) 



< 



vV + 1 - s 
s 2 + l 



) 
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We may assume that s' > 1, therefore 

/ Vs' 2 + 1 - s' \ / 1 \ m 

It is straightforward to verify that for s > 1 

\J s 2 + 1 - s > 

OS 

Therefore, since we can assume n > 1, it follows that 

v^TT-j ^ = _JL JL_ _ 1 

s 2 + 1 > s 2 + 1 ~ 3s 3 + 3s > 6s 3 ~ 6 • 10 6 n 6 

Hence, by choosing s' > 10 3 n 3 v / 6m, the second constraint is fulfilled. 
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